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ABSTRACT 

There are irreducible differences between the Hubble constant measured 
locally and the global value. They are due to density perturbations and finite 
sample volume (cosmic variance) and finite number of objects in the sample 
(sampling variance). We quantify these differences for a suite of COBEl- 
normalized CDM models that are consistent with the observed large-scale 
structure. For small samples of objects that only extend out to 10,000 km/sec, 
the variance can approach 4%. For the largest samples of Type la supernovae 
(SNela), which include about 40 objects and extend out to almost 40,000 
km/sec, the variance is 1 — 2% and is dominated by sampling variance. Sampling 
and cosmic variance may be an important consideration in comparing local 
determinations of the Hubble constant with precision determinations of the 
global value that will be made from high-resolution maps of CBR anisotropy. 

Subject headings: cosmology: distance scale, theories, observations 
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1. Introduction 

On the largest scales 100 Mpc) the Universe is well described by homogeneous and 
isotropic expansion satisfying the Priedmann equation. The global rate of the expansion at 
the current epoch is defined to be the Hubble constant (= Hq), the fundamental parameter 
of cosmology that sets the size and age of the Universe. Many methods have been used to 
measure Hq (see, e.g., Preedman 1996). Currently, the use of SNela as standard candles 
yields the smallest estimated measurement error (Riess, Press and Kirshner 1996; Hamuy 
et al 1996; Saha et al 1996). 

On small scales 100 Mpc), the Universe is significantly inhomogeneous. Because 
density fluctuations give rise to deviations from isotropic and homogeneous expansion 
(peculiar velocities) the expansion cannot be characterized by a universal expansion rate 
and measurements within a small, finite region will yield a local expansion rate (= Hi^) 
which is not identical to the global expansion rate. The difference arises from two factors: 
the finite physical size of the sample (Turner, Cen & Ostriker 1992; Nakamura and Suto 
1995; Shi, Widrow and Dursi 1996; Turner 1997) and the limited number of objects in 
the sample. Because of peculiar velocities, the average expansion rate in a finite volume 
is different from the global expansion (cosmic variance). Moreover, because only a small 
number of points within the volume are sampled the expansion rate defined by these points 
can deviate from the average expansion rate for the volume (sampling variance). These 
effects are different from measurement error, which can be reduced by better measurements 
and/or better standard candles. Cosmic and sampling variance can only be reduced by 
increasing the sample volume and the sampling density. 

In this Letter, we quantify cosmic and sampling variance for two samples of SNela and a 
cluster sample with TuUy-Pisher distances. In so doing we use a suite of COBE-normalized 

CDM models that are consistent with the measurements of the level of inhomogeneity on 
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scales less than about 300 Mpc and which should therefore provide a reasonable estimate 
for the variance that arises due to inhomogeneity in the Universe. 

2. Methodology 

The deviation of the local expansion rate measured by an observer at position r from 
the global Hubble constant is given by 



where v(r' — r) • (f — f) is the measured radial peculiar velocity at r' and W{r' — r) is 
the window function that characterizes the sample volume and sampled points within the 
volume (more below). The quantity v(r' — r) • (f — f) consists of two parts, the actual 
radial peculiar velocity and its measurement error, which is roughly the measurement error 
of the distance scaled by Hq. The measurement error can be shrunk by improving distance 
measurements, but the real radial peculiar velocity is an intrinsic deviation from the Hubble 
flow that is determined by the underlying density fluctuations. It is its contribution to 
5H{y)/Hq that is irreducible. 

The peculiar-velocity fleld depends upon the underlying power spectrum of density 
perturbations, which in turn, depends upon the cosmological scenario. We shall investigate 
a number of CDM models. In linear-perturbation theory. 



where 5{v') is the density fluctuation [p(r') — p]/ p and is the fraction of critical density 
in matter that clusters. Taking the Fourier transform we flnd. 
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where 5(k) = (27:)'^^'^ J d^r5{r)e^^ '' and Z(k) = J d^rW{r)r/r e^^ ''. The variance of 
5H{v)/Hq is 

For a top-hat, spherical window function, 

_ o sin(A:i^) -/o^^da;sin(x)/x 

where W{y' — r) = Q{R — |r' — r|)/(47ri?''^/3), Q{x) is the step function and R is the radius 
of the top-hat sphere (Shi, Widrow and Dursi 1996). Since the top-hat spherical window 
function samples every point within the sphere, {{5H{y) / HqY)^/'^ reflects only the cosmic 
variance associated with the finite volume of the spherical sample. 

Real data sets do not sample every point in space, rather they sample a number of 
objects with positions and redshifts Zq. The radial peculiar velocities of these objects are 

V5 • fg = CZq - HoVg (6) 

with uncertainties cr^ which are essentially the uncertainties in distances scaled by 
Hq because the uncertainties in Zg are relatively small. Random motions due to local 
nonlinearities, characterized by a one-dimensional standard velocity dispersion a*, may be 
added to cr^ in quadrature. Typically cr* ~ 10^ km/sec, and is therefore negligible when 
the sample depth is as large as ~ 10^ km/sec. Corrections to the linear Hubble law due to 
deceleration can also be safely ignored for samples with z <^ 1. With these approximations 
it follows that (Shi 1997) 



where 
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indices q, q' and q" denote summation over objects, and indices i, j, I, m are spatial indices 
that run from 1 to 3. Now {{5H{v) / HqYY/'^ includes both cosmic variance and sampling 
variance. 

3. Results 

Using Eqs. (|I|) to (|^) we have calculated the cosmic variance portion of {{5H{v) / H^YY/"^ 
with R =7,000, 10,000, 15,000, 20,000, 25,000, and 30,000 km/sec (see Fig. 1). The 
underlying cosmological models are a suite of COBE-normalized CDM models that are 
consistent with large-scale structure on scales from about 300 Mpc to 0.1 Mpc (Dodelson, 
Gates, and Turner 1996). The CDM models include a model with a low value of the Hubble 
constant, with significant tilt, with 20% light neutrinos, with additional radiation, and 
with a cosmological constant. In addition, we have included an open CDM model that 
is consistent with large-scale structure measurements and for completeness, a standard 
CDM model, which has excessive inhomogeneity on scales less than 300 Mpc. These models 
should serve well to span theoretical expectations for the level of inhomogeneity on the 
scales relevant for local variations in the Hubble constant. Their cosmic parameters are 
summarized in Table 1. 

The cosmic variance of at i? = 7, 000 km/sec is significant, ranging from about 2% 
to almost 4%. At i? = 10,000 km/sec it has fallen to 1% to 2%, and quickly declines to 
below 1% at i? = 15, 000 km/sec for all models. At a depth of i? = 30, 000 km/sec, which is 
reached by SNela, the cosmic variance is only about 0.2%. 

Using Eq. (H) we have calculated {{5H{y) / Hof'Y^'^ for volumes that are sampled by 
a finite number of points, so that both cosmic and sampling variance are included. For 
definiteness we use the SNela sample of Riess et al (1997), for which = 65 km/sec/Mpc, 
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the SNela sample of Hamuy et al (1996), for which //^ = 63.1 ± 3.4 ± 2.9 km/sec/Mpc, 
and the TuUy-Fisher sample of 36 clusters used for the template TuUy-Fisher relation in the 
Mark III catalogue (Willick et al 1997). Our results are compiled in Table 2. 

The intrinsic variance of iJ^ measured in the two SNela samples is around 1%, far less 
than measurement error. Although the SNela sample of Riess et al (1997) is deeper and 
has more objects than that of Hamuy et al (1996), its effective depth is not as large because 
many of the SNela are nearby. Due to its shallow depth, the TuUy-Fisher cluster sample 
has a larger intrinsic variance, between 1.5% and 3%. 

Figure 2 illustrates the effect of finite samphng for the SNela sample of Riess et al 
(1997). Comparing this plot to Fig. 1, it can be seen that the cosmic + sampling variance is 
more than twice the cosmic variance. At present, sampling variance dominates the intrinsic 
variance for the SNela samples. Since sampling variance scales roughly as the inverse square 
root of the number of objects, it can be shrunk to less 1% for all viable CDM models if the 
number of SNela is doubled. 

4. Summciry 

There are intrinsic and irreducible differences between the locally measured value of 
the Hubble constant and the global value. They arise due to finite sample volume (cosmic 
variance) and finite sample size (sampling variance) and can of course be of either sign. 
Cosmic variance and sampling variance cannot be reduced by better measurements or 
better standard candles. 

We have calculated the theoretical expectations for a suite of COBE-normalized CDM 
models that are consistent with measurements of large-scale structure on the scales that 
give rise to the cosmic variance portion. For samples that only extend out to 7,000 km/sec 
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the cosmic variance alone can be close to 4%; for samples of around 30 objects that extend 
out to 10,000 km/sec cosmic + sampling variance is between 2% and 4%. For samples of 
around 40 objects that extend out to 40,000 km/sec the the total variance is between 0.5% 
and 1.5%, with the cosmic variance contribution being less than 0.25%. 

As local measurements of the Hubble constant become more precise, cosmic variance 
and sampling variance will become a larger portion of the error budget and may be 
important when comparing local measurements with the better than 1% determinations of 
the Hubble constant anticipated from high-resolution maps of CBR anisotropy (Jungman 
et al 1996). 

The authors thank Adam Riess for providing the positions of their unpublished SNela. 
X.S. is supported by grants NASA NAG5-3062 and NSF PHY95-03384 at UCSD. M.S.T. is 
supported by DoE (at Chicago and Fermilab) and by the NASA through grant NAG 5-2788 
at Fermilab. 
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Table 1. Parameters of seven CDM models* 





sCDM 


h = 0.4 CDM 


tCDM 


i/CDM 


rCDM 


ACDM 


oCDM 


Otot 


1 


1 


1 


1 


1 


1 


0.4 


Om 


1 


1 


1 


1 


1 


0.4 


0.4 


Ob 


0.1 


0.16 


0.08 


0.07 


0.08 


0.06 


0.07 













0.2 











h 


0.5 


0.4 


0.55 


0.6 


0.55 


0.65 


0.6 


n 


1 


1 


0.7 


1 


0.95 


1 


1.1 



*Models and are from Dodelson, Gates and Turner (1996). Their power spectra are 
based upon Bardeen et al (1986) (transfer function), Bunn and White (1997) (COBE 
normahzation), Sugiyama (1995) (effect of fis on transfer function), Ma (1996) (z/CDM 
transfer function). The oCDM model is White and Silk (1996) (oCDM). With the 
exception of standard CDM, which is included only for completeness, all models are 
consistent with measures of large-scale structure on scales from 300 Mpc to 0.1 Mpc. 
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Table 2. Cosmic + sampling variance for three samples. 





SNela 


SNela 


TuUy-Fisher 


Reference 


Riess et al 1997 


Hamuy et al 1996 


Willick et al 1997 


Maximal depth 


37,000 km/sec 


30,000 km/sec 


11,000 km/sec 




~21,000 km/sec 


~ 27, 000 km/sec 


~ 10, 000 km/sec 


# of objects 


44 


26 


36 


sCDM 


1.4% 


1.4% 


3.1% 


h = 0.4 CDM 


1.1% 


1.1% 


2.6% 


rCDM 


1.0% 


0.9% 


2.3% 


uGBM 


1.3% 


1.3% 


3.1% 


tCDM 


0.8% 


0.8% 


1.8% 


ACDM 


1.0% 


0.9% 


2.4% 


oCDM 


0.7% 


0.6% 


1.6% 



*^peak is the wave number where |Z(k) • reaches maximum. 
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Figure Captions 

Fig. 1: Cosmic variation as a function of the radius of the sample volume for the seven 
COBE-normalized CDM models considered. 

Fig. 2: Cosmic + sampUng variance for the SNela sample of Riess et al (1997). The 
numbers indicate the number of objects within the spherical sample volume. 
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